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1  Introduction  and  a  reorientation 

Contract  F49620-92-J-0056  dealt  with  mathematical  aspects  of  “electromagnetic  bullets.” 
The  research  had  tw;  objectives: 

la)  Given  a  solenoidal  vector  field  J  defined  throughout  a  region  R,  find,  if  possible, 
a  set  of  Clebsch  Poicntis«-l3[27]  rj  such  that 

J  =  grad  4  X  grad  77  in  R  (1) 

Note:  If  v/e  pick  a  point  p  in  R  then  Clebsch  potentials  are  known  to  exist  in 
some  neighborhood  U  of  p. 

lb)  In  general,  find  the  topological  obstructions  to  finding  a  set  of  globally  defined 
Clebsch  potentials. 

2)  Investigate  the  (time  domain)  techniques  required  to  fonnxilate  a  nonlinear 
version  of  the  Lorentz  Reciprocity  Law. 

From  the  outset,  I  mention  that  much  folklore  appearing  in  the  published  literature  has 
serious  errors  which  must  be  articulated  in  order  to  achieve  the  objectives  of  the  proposed 
research.  Specifically, 

» 

1.  It  is  widely  believed  (Kuznetsov  and  Mikhailov[54]  are  often  dted  as  a  source)  that 
a  set  of  global  Clebsch  potentials  can  be  found  by  “pulling  back”  the  area  formed 
on  the  sphere  S^.  In  the  physics  literature  this  is  called  an  n-field  representation. 
Although  the  underlying  intuition  can  be  obtained  by  appealing  to  the  gyrovector 
in  micromagnetics [50,  51],  I  have  8hown[53]  that  the  gyrovector  construction  is  not 
possible  for  a  general  solenoidal  vector  field.  (The  present  research  follows  up  on  this.) 

2.  The  nonlinear  version  of  Lorentz  Reciprocity  cannot  be  articulated  in  terms  of  fre¬ 
quency  domain  concepts.  Similarly,  the  fact  that  a  Fourier  Transform  relates  an  an- 
teuna’s  far  field  to  its  sources  cannot  be  used  to  “eicplain”  why  electromagnetic  bullets 
or  missiles  cannot  exist  since,  by  construction,  there  is  no  far  field  for  a  bullet.  Over 
the  last  year  it  has  become  clear  that  one  has  to  deal  with  Maxwell’s  Equations  as  a 
system  of  hyperbolic  p.d.e.’s  and  avoid  the  temptation  of  using  elhptic  theory  which 
is  applicable  when  taking  a  Fourier  Transform  (as  engineers  are  trained  to  do)  and 
playing  with  Helmholtz’s  equation.  The  way  to  achieve  these  goals  is  to  reexamine  the 
“raison  d’etre”  for  the  use  of  the  Radon  transform  in  these  hyperbolic  problems. 

In  reference  to  this  second  point,  I  have  discovered  that  the  Mathematics  Literature  has  a 
natural  reason  for  the  appearance  of  the  Radon  transform[74,  5].  Specifically,  if  one  tries 
to  describe  the  lacunas  (Latin  for  hole)  of  fundamental  solutions  of  hyperbolic  p.d.e.’s,  one 
notices  that  the  support  (in  space-time)  of  the  singulu  part  of  a  solution  lies  on  a  light  cone 
made  of  rays.  These  rays  are  identified  with  the  points  of  a  projective  space  which  form  the 
domain  of  integration  in  the  Radon  transform.  Although  not  part  of  the  orifpnal  research 
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plan,  pursuit  of  this  observation  seems  crucial  for  a  deeper  understanding  of  the  bullet  idea 
and  of  reciprocity. 

In  addition  to  the  ties  with  lacunas  of  hyperbolic  equations,  my  research  has  focused 
primarily  on  two  areas: 

1.  The  role  of  higher  order  linking  of  flux  lines  as  an  obstruction  to  Clebsch  potentials, 
(i.e.  Massey  products  in  a  suitable  cohomology  ring  come  to  play.) 

2.  Finding  an  underlying  incredible  Hauniltonian  system  with  enough  structure  to  indicate 
obstructions  to  globally  defined  Clebsch  potentials.  (So  far  we  have  three  such  systems 
but  the  clearly  articulated  obstructions  are  still  forthcoming.) 

The  remainder  of  the  report  develops  the  two  above  topics  with  emphasis  on  the  first. 

2  Massey  Products  and  obstructions  to  Clebsch  po¬ 
tentials 

2.1  Relationship  to  recent  developments  in  algebraic  topology 

Massey  Products  [56,  57,  80]  are  a  way  of  articulating  higher  order  linking.  For  example,  the 
curves 


are  unseparable  because  they  are  linked.  Higher  order  linking  can  be  illustrated  by  Bor- 


a  previous  publication[48]  I  have  stumbled  across  a  need  for  Massey  products  in  another 
context.  This  higher  order  linking  is  important  in  the  context  of  Clebsch  potentials  since 
they  describe  an  obstruction  to  having  a  set  of  single  valued  globally  defined  potentials.  To 
develop  this  connection  further,  recall  that  in  a  neighborhood  U  of  a  point  p,  a  solenoidal 
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yector  field  J  is  expressible  in  terms  of  Clebsch  potentials  and  hence,  it’s  vector  potential 
(which  we  denote  by  F)  has  a  Monge  potential  representation: 

J  =  curl  F  where 

F  =  ^  grad  +  grad  inU  ' 

Ties  between  Clebsch  potentials  and  Massey  products  are  best  motivated  by  appealing  to 
the  “helicity  of  the  vector  potential.”  The  helidty  of  a  vector  field  G  is  ^ven  by 

y  G  •  curl  G  dV  (3) 

and  its  interpretation  in  terms  of  twisted  and  tangled  flux  lines  has  been  developed  by  Moffatt 
in  the  context  of  magnetohydrodynamics[63,  64,  65,  66,  67,  68,  69,  70,  71,  72,  2].  Berger 
and  Field[9]  have  continued  this  investigation  and  Berger  has  given  explicit  expressions  for 
Massey  Products  in  the  notation  of  pedestrian  vector  analysis  [10].  In  my  publication[53]  I 
show  that  an  obstruction  to  globally  defined  Clebsch  potentials  follows  from  the  fact  that 
the  helicity  of  any  vector  potential  described  by  global  Monge  potentials  vanishes.  More 
concretely,  nonzero  total  helicity  implies  an  obstruction. 

Further  insight  into  the  tie  between  helicity  and  the  global  obstructions  to  Clebsch  po¬ 
tentials  can  be  found  by  appealing  to  the  Frobenius  integrability  condition: 

If  F  •  curl  F  =  0  in  U 

then  F  =  ^  grad  tf  for  some  ^  and  rj 

» 

A  quick  calculation  shows  that  the  converse  of  the  implication  (4)  is  true.  By  the  gauge 
freedom  of  V')  apparent  in  (2),  we  see  that  setting  ^  equal  to  zero  ensures  that  the  helicity  of 
F  can  be  made  to  vanish  in  some  neighborhood  of  any  point  p!  Since  the  helicity  integral  (3) 
subject  to  boundary  constraints  is  usually  a  fixed  nonzero  number  we  see  that  the  obstruction 
to  globally  defined  Monge  (and  hence  Clebsch)  potentials  follows  from  the  twisting  and 
tangling  of  flux  lines  and  from  Massey  products  describing  “linking”  of  closed  flux  lines. 

The  tie  between  helicity  and  Massey  products  or  other  obstructions  to  globally  defined 
Clebsch  potentials  is  still  incomplete  but  recent  mathematical  developments  point  to  a  deep 
connection  between  these  ideas.  To  see  this  we  note  that  all  of  the  above  arguments  can  be 
rephrased  in  terms  of  differential  forms  in  order  to  demonstrate  that  they  are  independent 
of  the  Riemannian  metric  on  the  space.  This  coordinate  free  approach  then  ties  into  some 
basic  facts  and  some  recent  history: 

1.  There  is  a  non-abelian  generalization  of  the  helicity:  the  “Chem-Simmons  Secondary 
Characteristic  Classes,”[26, 73].  That  is,  if  one  thinks  of  a  1-form  as  a  U(l)  connection 
and  it’s  exterior  derivative  as  a  curvature,  the  helicity  reduces  to  a  Chern-Simmons 
form.  More  generally,  if  .4.  is  a  Lie  algebra  valued  1-form,  and 

dA-\-A/\A  (5) 

the  associated  curvature,  the  Chem-Simmons  form  looks  like: 

j  ir[A^dA■\■'^A^A^A]  (6) 


4 


Ed  Witten[93]  has  constructed  an  associated  topological  quantum  field  theory  which 
yielded  the  Jones  polynomial  of  knot  theory  [44]  and  an  extension  to  three-mamfolds 
(as  well  as  a  fields  medal  for  Witten).  The  underlying  philosophy  is  that  if  one  starts 
with  a  metric  invariant  Lagranian  (the  Chem  Sinunons  form  in  this  case)  then  the 
expectation  values  associated  with  any  metric  independent  observable  are  necessarily 
topological  invariants. 

2.  To  most  mathematicians,  Witten’s  use  of  Feynmann  path  integrals  is  tantamount  to 
witchcraft,  and  so  an  independent  approach  to  the  Jones  polynomial  is  sought  via 
orthodox  algebraic  topology.  Such  an  understanding  is  beginning  to  emerge  through 
the  work  of  Vassilliev[89,  12]  who  uses  the  Wderstrauss  approximation  theorem  to 
build  a  sequence  of  spaces  approximating,  and  in  some  sense  converging,  to  the  “Space 
of  embeddings  of  S^  into  (i.e.  the  space  of  aU  knots).  Vassiliev’s  invariants  are 
then  extracted  from  a  spectral  sequence  involving  the  reduced  homology  groups  of 
these  spaces.  Birman  and  Lin  have  extracted  the  Jones  polynomial  from  Vassiliev’s 
invariants.  See  Birman[12]  for  a  recent  survey  article  (especially  Section  Four).  Finally, 
it  is  important  to  say  that  the  entire  formalism  goes  over  to  the  case  of  links  (i.e. 
embeddings  of  several  copies  of  S^  in  M^). 

3.  What  is  most  exciting  is  the  emer^g  relationship  between  perturbative  approaches 
to  the  Chern-Simmons  topological  quantum  field  theory  and  Vassiliev’s  invariants. 
The  recent  work  of  Baez[6]  and  Kontsevich[47],  appear  to  be  the  first  successes  at 
making  the  Feynmann  integral  approach  to  Vassiliev  invariants  a  psurt  of  orthodox 
algebraic  topology.  The  bridge  to  link  invariants  such  as  Massey  products  is  also 
emerging  through  attempts  to  decompose  Vassiliev’s  invariants  into  a  sum  containing 
well  known  invariants  and  “other  terms.”  The  Milnor  ^-invariants  are  known  to  arise 
in  this  way [7,  13].  In  the  context  of  knot  complements,  earlier  work[75,  88]  has  shown 
that  these  /i-invariants  are  just  Massey  products. 

4.  Another  framework  for  Jones-like  invariants  comes  from  the  study  of  “exactly  solvable 
problems  in  statistical  mechanics,”  R-matrices,  and  quantum  groups.  Although  it  is 
beyond  my  competence  to  describe  these  developments,  it  is  important  to  point  out  one 
connection  with  my  previous  work  on  the  helidty  functional.  In[49]  I  consider  the  fi¬ 
nite  element  discretization  of  the  helidty  functional  (3)  by  means  of  Whitney  forms  [90]. 
Over  a  tetrahedron  A,  the  discretized  expression  involves  six  independent  edge  vari¬ 
ables  hy,  1  <  i,  j  <  4,  where  the  indices  correspond  to  vertices  of  the  tetrahedron  and 
hy  =  -hjj.  Upon  discretization  we  have 

f  H  •  curl  H  dV  — *  hi2h34  —  hi3h24  +  hi4h23  (7) 

/A 

The  right  hand  side  is  the  PfafRan  of  the  skew-symmetric  4  by  4  matrix  formed  by 
the  hij  and  it  clearly  preserves  the  metric  independence  of  helidty.  This  expression 
has  the  same  symmetries  as  the  **6-j  symbols"  which  appear  in  the  tensor  products  of 
group  representations  and  hence  seem  to  point  to  a  more  direct  Knk  with  guantum  group 
invariants.  Although  intriguing,  a  more  articulate  description  of  what  is  going  on  is 
not  apparent  to  me  or  prominent  experts  whom  I  have  approached. 
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The  above  four  points  are  summarized  schematically  in  Figure  1.  Our  objective  is  to 
articulate  how  “linking  and  tangling  of  flux  lines”  provide  a  global  obstruction  to  defining  a 
set  of  global  Clebsch  potentials.  Noting  that  Monge  potentials  can  always  be  chosen  so  that 
the  helicity  density  vanishes  locally,  the  helicity  integral  provides  an  obstruction  to  globally 
defined  Clebsch  potentials  and  a  bridge  to  an  algebraic  description  of  the  obstruction  (e.g. 
Massey  products). 


Figure  1:  Question:  Can  Whitney  forms  and  finite  elements  provide  a  more  direct  link? 
Possible  Answer:  Pfaifian  Expression  for  Helicity 


2.2  A  diversion  on  Massey  products  and  Morse  inequalities 

It  is  my  belief  that  Massey  products  have  a  fundamental  role  to  play  in  critical  point  theory 
via  the  “Mountain  Pass  Lemma”  and  its  generalizations  [81].  The  basic  underlying  intuition 
is  illustrated  by  issues  in  biomedical  imaging[45,  78].  Formally,  we  would  like  to  make  a 
connection  between  Morse  theory  and  Massey  products.  Let  us  be^n  by  reviewing  the 
fundamentals  of  Morse  theory[61,  14,  34,  62]. 

Consider  a  compact  n-dimensional  manifold  M,  denote  the  i^^  Betti  number  (i.e.  the  rank 
of  the  ith  homology  group)  by  A(M),  and  form  the  Poincare  polynomial: 

(8) 

k=0 

The  critical  points,  {p,},  of  a  function  f,  mapping  M  to  are  said  to  be  nondegenerate 
if  the  Hessian  matrix  (i.e.  matrix  of  second  partial  derivatives)  has  full  rank  at  every  critical 
point.  A  function  with  nondegenerate  critical  points  is  called  a  Morse  function.  (Morse 
functions  are  not  rare  in  that  they  form  an  open  dense  set  in  the  space  of  all  functions.) 
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Given  a  Morse  function  on  M,  we  define  the  index  of  the  critical  point,  pj,  to  be  the  number 
of  negative  eigenvalues  of  the  Hessian  matrix  evaluated  at  pj.  (Note:  the  eigenvalues  are  real 
because  the  Hessian  is  symmetric,  and  nonzero  since  pj  is  nondegenerate.)  Let  mi  be  the 
number  of  critical  points  of  index  i,  and  form  the  Morse  polynomial: 

M.{{)  = 'ti'‘ mt({)  (9) 

k=0 

A  fundamental  result  of  Morse  theory  asserts  that 

M(f)-Pt(M)  =  (l  +  t)Qt(f)  (10) 

where  Q(f)  is  a  polynomial  with  nonnegative  integer  coefficients.  An  equivalent  formulation 
(the  Morse  inequalities)  is  found  by  equating  powers  of  t,  solving  for  the  coefficients  of  Qt(0 
SLud  using  the  iionnegativity  of  these  coeffidents  to  obtain  the  inequalities: 


mo 
mi  —  mo 
mj  —  mi  -|-  mo 


>  /3o 

>  01-00 

^  02  —  01+00 


:  > 


(11) 


Given  the  Morse  Inequalities  expressed  in  the  form 'of  equation  (10),  two  strategies  for  their 
use  emerges: 


1.  If  Pt(M)  is  known,  then  one  has  a  lower  bound  of  the  number  of  critical  points  of  every 
index.  Spedfically,  a  very  loose  bound  is  given  by  m;  >  0i  for  all  i. 

2.  If  no  consecutive  powers  of  t  occur  in  the  Morse  polynomial  then  Qt(I)  =  0  and  so 
Alt(f)  =  Pt(M).  In  this  case  knowledge  of  the  critical  points  gives  complete  knowledge 
of  the  Betti  numbers. 


A  function  for  which  Mt{i)  is  equal  to  Pt(M)  is  called  a  perfect  Morse  function.  In 
general  there  are  topological  obstructions  to  finding  a  perfect  Morse  function.  For  example, 
the  Morse  inequalities  are  also  valid  if  homology  groups  are  computed  using  coeffidents  from 
any  field  (e.g.  where  p  is  a  prime).  Hence  a  perfect  Morse  function  is  only  possible  when 
the  homology  groups  of  M  are  torsion  free. 

With  this  background  we  can  return  to  the  topic  of  Massey  products.  Concretely  I 
would  like  to  conjecture  that  the  obstruction  to  finding  a  perfect  Morse  function  depends  on 
Massey  products  (just  as  the  obstruction  to  globally  defined  Clebsch  potentials  seems  to). 
Specifically,  I  have  a  conjecture: 

Q  _  Q  ^  /  All  Massey  products  in  the 
~  (  cohomology  ring  of  M  are  trivial 

I  have  bits  of  circumstantial  evidence  to  support  this  conjecture.  I  will  present  them  before 
I  proceed  to  outline  how  I  may  go  about  proving  (12). 
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1.  The  construction  of  perfect  Morse  functions  plays  a  great  role  in  computing  Poincare 
pol3moinials  in  various  contexts.  For  example  Yang-Mills  Theory  on  Riemann  Surfaces, 
stable  homotopy  groups  of  Lie  groups[14],  Grassmannians,  etc.  These  examples  seem 
to  indicate  that  the  above  conjecture  is  overly  ambitious.  However,  in  these  cases  the 
underlying  manifold  has  a  Kahler  structure,  and  it  is  known  that  Kahler  manifolds 
have  vanishing  Massey  Products  [39].  Hence,  in  order  to  make  progress  with  the  above 
copjecture,  we  can  search  the  literature  for  answers  to  the  following  two  questions: 

(a)  Do  all  Kahler  manifolds  admit  a  perfect  Morse  function?  If  not  what  are  the 
obstructions? 

(b)  Are  there  classes  of  non-Kahler  manifolds  which  are  known  to  admit  perfect  Morse 
functions? 

2.  In  both  the  calculus  of  variations  and  finite  dimensional  critical  point  theory,  the 
“mountain  pass  lemma”  is  frequently  used  to  exhibit  saddle  point  extrema.  (This  is 
a  powerful  technique  in  the  theory  of  nonlinear  elliptic  equationsfSl]).  Underlying 
the  proof  of  the  mountain  pass  lemma  is  a  linking  argument  which  leads  to  deeper 
results.  The  notion  of  linking  describes  the  simplest  of  cohomology  operations  and 
Massey  products  appear  to  be  the  right  framework  for  generalized  results  in  critical 
point  theory. 

To  see  how  Massey  products  can  come  into  Mgrse  Theory  naturally,  we  have  to  recall 
how  the  polynomial  Qt(f)  arises  in  the  first  place.  The  level  sets  of  a  function  f:  M  -+  jR 
define  a  “stratification”  of  the  manifold  M.  Let  us  perturb  f  such  that  all  critical  points  are 
nondegenerate  and  each  level  set  contains  at  most  one  critical  point.  Furthermore,  define 
submanifolds 


Xc  =  {x  G  M  I  f(x)  <  c}  (13) 

If  denotes  the  A-dimensional  unit  ball  whose  boundary  is  the  (A  —  l)-dimensional  sphere 
and  'v/  is  used  to  denote  a  homotopy  equivalence,  then  it  is  a  basic  result[14,  34,  61], 
that  for  sufficiently  small  e  >  0,  we  have  the  following  “change  in  topology”  formula: 


Xc+e- 


Xc-e 
Xc-*  U 


if  c  is  not  a  critical  value 

if  c  is  a  critical  value  tied  to  a  critical  point  p 

of  index  A. 


(14) 


To  see  how  the  polynomial  Qt(f)  emerges  we  now  appeal  to  algebraic  topology.  Consider 
the  above  “change  in  topology  formula”  for  each  critical  point  pj  and  corresponding  critical 
value  Cj.  Introduce  the  inclusion  and  projection  maps: 


ij  :  Xc^._e 
ITj  : 


Xc3+e/Xc3-e 


(15) 


and  consider  the  long  exact  sequence  of  homology  groups  arising  at  the  j^  critical  point: 
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where  ij,i  and  TTjj  are  induced  mappings  on  the  1th  homology  groups  and  is  the  “connecting 
homomorphism.”  If  we  define  the  Poincare  polynomial  of  the  connecting  homomorphism  as 

P(Im  5j)  =  dim(Im  t"^  (17) 

k=0 

then,  by  standard  techniques [34],  we  have 

=  P{X,,+.)  -  P(X.,..)  +  (1  + 1)  P(Im  «i)  (18) 

We  win  now  outline  how  the  Morse  inequalities  follow  if  we  extract  the  algebraic  content 
of  equation  (14)  and  consider  equation  (18)  for  each  critical  point  and  sum.  First,  from 
equation  (14)  we  deduce  a  homotopy  equivalence 

Xcg_,j+e  ~  Xcj-e  (^9) 

and  so,  by  the  homotopy  invariance  of  homology  groups,  we  have 

Hk(Xcy_„+.)  Hu(Xc,-e)  for  all  k 
and  Pt(Xc(,_.,+e)  =  Pt(Xc,-.) 

Notice  that  if  f  has  N  critical  points  Xcq-c  is  the  empty  set  and  X^^^-e  is  M.  Equation  (20) 
then  shows  that  if  we  consider  a  copy  of  equation  (18)  for  each  critical  point  and  sum,  the 
first  two  terms  on  the  right  hand  side  form  a  telescoping  sum  so  that 

|:Pt(Xc,+e,Xc,-e)  =  Pt(M)  +  (1  +  t)2Pt(Im  Si)  (21) 

j=l  j=l 

Next,  we  continue  to  extract  algebraic  consequences  of  equation  (14).  As  we  cross  a  critical 
value,  in  equation  (14)  is  attached  to  Xq-e  along  an  and  this  provides  us  with  a 
relative  homotopy  equivalence  which  descends  to  an  isomorphism  of  homology  groups  and 
equality  of  Poincare  polynomials: 

(Xq+e,Xcj_e)  ~  (S^‘,pt.)  implies  1 

Hic(Xc,+«,Xcj-.)  Hk(S^,pt.)forallk,  \  (22) 

and  so  Pt(Xq+«,Xcj_,)  =  Pt(S^,pt.)  =  t^  J 
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4. 


Substituting  this  expression  for  the  Poincare  polynomial  into  the  left  hand  side  of  equa¬ 
tion  (21)  and  collecting  terms  with  same  index,  the  left  hand  side  is  seen  to  be  the  Morse 
polynomial  defined  by  equation  (9).  Hence 

Mt(f)  =  P(M)  +  {1  + 1)  |:P(Im  dj)  (23) 

j=i 

At  this  point  we  recover  the  Morse  inequalities  (10)  and  Qt(f)  can  be  identified  with  the 
summation  on  the  right  hand  side  -  but  we  can  do  better! 

Given  the  homology  groups  of  the  sphere  modulo  a  point,  equation  (22)  shows  that 
Hk(Xcj+*,Xcj_e)  is  trivial  unless  k=Aj.  Hence,  in  the  exact  sequence  (16),  the  connecting 
homomorphisms  are  all  trivial  except  for  possibly  This  means  equation  (17)  can  be 
simplified  to: 

Pt(Im  tfj)  =  dim(Im  ^j.Aj)t^^"^  (24) 

Substituting  this  into  (23)  and  identifying  Qt(f)  ^ves 


(25) 


At  this  point  we  can  return  to  the  conjecture  (12).  The  conjecture  amounts  to  saying  that  if 
we  compute  the  cohomology  of  the  manifold  M  by  using  a  cell  decomposition  derived  from 
the  critical  points  and  level  sets  of  a  Morse  function,  then  nontrivial  Massey  products  in 
the  cohomology  ring  should  force  dim(Im  5j,Aj)  to  be  nonzero  for  some  collection  of  j’s.  To 
make  this  precise,  we  propose  to  compute  the  cohomology  of  M  by  appealing  to  a  spectral 
sequence  arising  from  the  cell  decomposition  of  M  and  to  reexpress  Qt(f)  in  terms  of  the 
differentials  of  the  spectral  sequence.  It  is  known  that  Massey  products  are  naturally  related 
to  the  differentials  of  a  spectral  sequence[59]  and,  hopefully  we  can  obtain  deeper  insight 
into  conjecture  (12)  by  pursuing  the  spectral  sequence  framework. 

Let  us  consider  two  other  contexts  which  shed  light  on  conjecture  (12): 


Q.(f)  =  Ef=.  dim(Im 


3.  Further  insight  into  conjecture  (12)  can  be  obtained  by  considering  the  obstruction  to 
finding  a  perfect  Morse  function  on  a  compact,  connected,  orientable,  three  dimensional 
manifold  M^.  To  start,  consider  two  embeddings.  Mi  and  M2,  into  JR^,  of  a  compact, 
connected,  two-dimensional  manifold  M,  of  genus  g.  The  embeddings  force  M,  (i  =  1, 
2)  to  be  orientable  and  hence  they  have  an  “inside”  and  an  “outside.”  Call  the  insides 
Hi  and  H2  so  that 

aHi  =  Mi  (i  =  l,2)  (26) 

The  Hi’s  are  called  handle  bodies.  If  we  have  a  dir somorphism  a  :  Mi  — »  Mi  then  we 
can  glue  Hi  to  Hi  along  their  boundaries  to  get  a  ihree  manifold  M^: 

M^  ~  Hi  U  Hi  (27) 
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It  is  a  basic  result  in  three  manifold  theory  that  any  connected,  compact,  orientable 
3-Manifold  admits  such  a  splitting  for  some  g[76].  This  is  called  a  Hegaard  splitting 
of  M^.  The  minimal  genus  g  for  which  a  ^ven  manifold  has  a  Hegaard  splitting 
is  a  topological  invariant  called  the  genus  of  and  is  denoted  here  by  g*(M^).  The 
genus  of  has  a  definition  in  terms  of  the  fundamental  group  of  M^. 

We  next  consider  Hegaard  splittings  which  arise  from  Morse  functions.  Consider  first 
the  following 

Basic  Fact:  Th“11.3(34l 


Let  M°  be  a  smooth  compact  connected  n-dimensional  manifold.  Then,  on 
M",  there  is  an  admissible  Morse  function  f,  where  0  <  f  <  n,  such  that  it 
has  only  cue  minimum  (point  of  index  0),  only  one  maximum  (point  of  index 
n),  and  each  critical  point  of  index  A  is  situated  on  the  surface  {~^(A). 

Note  that  this  statement  implicity  defines  the  notion  of  an  admissible  Morse  function. 
As  a  corollary,  we  see  that  an  admissible  Morse  function  defines  a  Hegaard  splitting. 
Concretely,  the  Morse  polynomial  of  an  admissible  Morse  function  f  :  M^  — » JR  looks 
like: 


M(f)  =  l  +  (l+t)tg(f)-J-t*  (28) 

where  g(f)  is  the  genus  of  the  two  dimensional  manifold  f~^(|)  and  the  “handle  bodies’* 
of  the  Hegaard  splitting  are  given  by: 


Hi  =  {x€M|f(x)<|} 

H2  =  {x  €  M  I  f(x)  >  |}. 

In  this  case,  the  Poincare  Polynomial  looks  like 

Pt(M'‘)  =  H-(l-ht)ti0i(M3)  +  t3 
Equations  (10),  (28),  and  (30)  enable  us  to  compute  Qt(f): 

(l  +  t)Qt(f)  =  Mt(f)-Pt(M=»)  =  (l  +  t)t(g(f)-A(M^)) 

i.e.  Q.(f)  =  t  (8(f)  -  ft(M^)) 

Equations  (25)  and  (31)  give  us  the  nontrivial  result: 


(30) 


(31) 


8(f)-A(M’)=  X;  dim(Im3u,)  (32) 
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By  definition,  for  any  admissible  Morse  function  we  have 


g(f)  -  A(M^)  >  g*(M=>)  -  >  0  (33) 

The  fact  that  g*(M^)  >  /9i(M^)  for  most  3-manifolds  shows  us  that  perfect  admissible 
Morse  fimctions  can  rarely  be  expected  to  exist.  Hence,  to  shed  light  on  conjecture 
(12)  we  can  consider  three  questions: 

•  Is  there  any  obstruction  to  finding  an  admissible  Morse  function  for  which 

g(l)  =  g-(M^)? 

•  Is  g’(M^)  —  i9i(M^)  expressible  in  terms  of  the  cohomology  ring  of  M^  and  in 
particular,  the  ‘^number”  of  nonzero  Massey  products?  (Valuable  insight  might 
be  obtained  by  relating  g*(M^)  to  Wi(M^)  to  Milnor’s  /^-invariants,  to  Massey 
products) 

•  How  can  alternate  definitions  of  g*(M^)  along  with  equations  (32)  and  (33)  be 
exploited? 

4.  Finally,  in  order  to  relate  Massey  products  to  Morse  theory,  we  may  try  to  relate 
Witten’s  proof  of  the  Morse  inequalities[92],  the  heat  equation  method  for  harmonic 
forms[30],  in  the  context  of  Witten’s  deformed  Lapladan,  and  John  Lott’s  work  on  the 
asymptotics  of  the  heat  kernel  and  its  relation  to  Massey  products[55]. 

This  connection  may  also  have  implications  for  molecular  biology [91,  85,  32]. 


3  Focus  of  ongoing  research 

3.1  Resolvable  intermediate  problems  regarding  Clebsch  poten¬ 
tials 

Surveying  the  literature,  it  is  clear  that  there  is  no  '^usable  answer”  to  the  question  of  finding 
the  topological  obstruction  to  having  globally  defined  Clebsch  potentials.  Furthermore,  the 
engineering  and  physics  literature  is  riddled  with  papers  which  are  completely  erroneous  when 
it  comes  to  articulating  the  global  existence  of  Clebsch  potentials.  On  the  other  hand,  the 
mathematics  literature,  in  particular  the  theory  of  G-8tructure8[25],  contact  structuresfS], 
and  foliations  contains  the  formalism  to  clearly  articulate  the  problem,  and  possibly  find 
a  solution  (even  if  it  is  crouched  in  algebraic  structures  which  are  not  easily  reduced  to 
numerical  computation).  For  this  resison,  I  have  formulated  several  intermediate  problems 
which  bridge  the  gap  between  the  classical  local  results  and  the  problems  addressed  by  the 
modem  formalism. 

Given  a  ‘‘nice”  region  R  in  and  a  divergence  zero  vector  field  G  in  R,  assume  that  we 
can  write 


G  =  curl  F  in  R 


(34) 


12 


i 


The  obstruction  to  doing  so  depends  on  the  cohomology  group  H^(R;JR)  which  we  assume 
to  be  trivial.  Next,  for  simplicity,  suppose 

F  ‘  curl  F  /  0  in  R  (35) 

We  can  then  describe  R  by  an  atlas  A  specified  by  a  Monge  potential  representation.  That 
is,  the  atlas  .4  is  a  collection  of  charts: 


l<i<n[A)  (36) 

where,  n(^)  is  the  number  of  charts  in  the  atlas  and,  by  the  definition  of  an  atlas,  the 
bijections  f\  satisfy 


(r.)-':  Ui-iR3 


Furthermore,  by  the  definition  of  Monge  potentials,  we  have 


(37) 


Taking  the  curl  we  have 


F  =  grad  1/5  +  grad  <jn  on  Ui 


(38) 


G  =  curl  F  =  grad  x  grad  rn  on  Ui  (39) 

and,  by  (35)  we  have: 

» 

F  •  curl  F  =  ^  0  on  Ui  (40) 

a(x,y,z) 

which  coincides  with  what  we  expect  from  the  definition  of  a  chart.  At  this  point,  it  is  very 
fruitful  to  pause  and  consider  the  geometric  consequences  of  what  we  have  done.  Equation 
(40)  implies 


while  equation  (6)  implies 


^(^i»^7i>^») 


1  on  Ui  n  Uj 


1  on  Ui  n  Uj 


(41) 


(42) 


From  this  vantage  we  see  that  there  are  several  interesting  and  relatively  simple  questions 
which  we  may  aisk. 


Question  1  Given  (R,  F)  as  above,  what  is  the  minimum,  over  all  possible  atlases, 
of  the  number  of  charts  in  an  atlas. 


Clearly  the  above  number  gives  a  measure  of  '*how  impossible”  it  is  to  find  a  set  of  global 
Clebsch  potentials  for  F  on  R. 

Equations  (41)  and  (42)  above  tie  our  problem  to  the  theory  of  G-structures,  and  the 
integer  which  answers  question  1  is  the  analogue  of  the  ‘‘Liustemik-Snirelman  Category”  for 
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some  cohomology  theory  cooked  up  from  the  constraints  on  the  differential  structure  on  our 
atlas. 

A  natural  variant  of  the  above  problem  is  the  following.  Consider  the  minimum  of  the 
answer  to  question  1  over  all  divergence  free  vector  fields  F  which  satisfy  equation  (35).  If  R  is 
a  closed  three-dimensional  manifold  then  this  yields  an  interesting  topological  invariant.  This 
invariant  is  intimately  related  to  the  work  of  Amold[2],  Moffatt[70,  71]  and  Freedman[36]  on 
using  helicity,  (mean  asymptotic  linking  number)  and  mean  asymptotic  crossing  numbers  to 
get  lower  bounds  on  the  energy  of  ideal  magnetohydrodynamic  flows. 

Another  intermediate  question  is: 

Question  2  Can  we  use  equation  (42)  to  try  to  build  a  differential  structure 

for  a  surface?  (Obstructions  would  also  be  obstructions  to  having  a  set 
of  globally  defined  Clebsch  potentials). 

If  so,  can  we  then  use  the  coordinate”  to  build  a  three 
dimensional  manifold?  That  is,  can  we  build  a  foliation  out  of 
surfaces  “0  =  const?” 

This  question  ties  back  to  the  use  of  Helicity  to  define  an  obstruction  to  globally  defined 
Clebsch  potentials. 

My  previous  work  in  micromagnetics[50,  51],  coupled  with  the  work  of  Amold[2,  4,  23, 
24,  46,  86],  Freedman[35,  36,  37],  Moffatt[65]  and  others,  indicate  that  a  fruitful  question  to 
consider  is  the  time  evolution  of  a  set  of  Clebsch  potentials.  Specifically,  in  the  ‘‘accessibility 
problem”,  where  a  vector  field  F  is  given  throughout  R  at  two  time  instants  ti  and  ta  and 
on  5R  for  t  €  [tj,  tj],  one  asks  if  there  exists  a  set  of  Clebsch  potentials  on  R  x  [ti,  t2].  The 
starting  point  of  such  an  investigation  is  to  derive  a  formula  which  says 


j  F  •  curl  F  dV  -H  (something)dS  dt  =  number  of 

“flux  lines  recoimected” 
in  R  X  [ti,t2] 

In  other  words,  let  I  =  [t^,  t2],  construct  some  3-form  u;  out  of  the  given  data  on  d(R  x  I), 
and,  without  the  hypotheses  given  by  equation  (35),  the  above  formula  says: 

/  w  =  oriented  flux  line  reconnections  in  R  x  I 
Ja(Rxi)  ^ 

3.2  Algebraic  gadgets  which  reveal  3-D  complexity 

The  obstruction  to  finding  a  set  of  global  Clebsch  potentials  is  ultimately  related  to  the 
twisting  and  linking  of  flux  lines,  and  the  task  at  hand  is  simplified  if  we  can  find  algebraic 
structures  which  articulate  this  behavior.  So  far  we  have  considered  linking  numbers  and 
Massey  products  in  cohomology.  I  have  found  that  a  useful  question  to  ask  is:  “How  un- 
contractible”  is  a  given  region  of  space?  It  turns  out  that  Massey  products,  K.T.  Chen’s 
iterated  integrals  [21],  loop  space  homology,  Denms  Sullivan’s  minimal  modeU[84],  the  lower 
central  series  of  link  groups,  and  Thurston’s  norm  on  homology[87]  all  come  to  play  and  are 
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intimately  related  in  this  context.  To  point  out  these  xdationships,  we  consider  the  Poincare 
homomorphism 


P  :  irx(R)  -  Hi(R;^)  (43) 

which  sends  the  homotopy  class  of  a  loop  in  R  to  its  homology  class.  The  map  is  onto  and 
its  kernel  is  the  commutator  subgroup  [xi,  xi].  Hence  we  have  the  isomorphism 


*•1  /  ki,xi]  ~  Hi  (44) 

This  isomorphism  plays  a  key  role  in  understanding  counterintuitive  aspects  of  making  “cuts” 
for  magnetic  scalar  potentials  in  three  dimensions  [40].  Our  immediate  concern,  is  to  find 
commutative  algebraic  gadgets,  more  complicated  than  Hi,  which  help  articulate  om- 
plexity  of  Xi.  (In  general,  for  xi’s  coming  from  three  dimensional  manifolds,  it  is  i  en 
known  whether  the  basic  group  theoretic  questions  ate  “decidable”  in  the  computer  &  icnce 
sense  of  the  word.)  We  begin  by  considering  the  graded  group  gr(xi)  constructed  out  of  xi, 
and  as  a  preliminary  step,  consider  the  lower  central  series  of  a  group  G,  defined  recursively 
by 


Gi  =  G,  Gn+i  =  (Gn,  G),  n  =  l,2, ...  (45) 

where  (•  ,  •)  denotes  a  commutator.  The  associated  graded  group  is  formally  defined  as 

gr(G)  =  ©“1  G'V  G'+'  (46) 

and  for  a  given  t  the  object  (see[41]) 

gr(G)  /  G‘+*  =  ©J^i  G'  /  G'+»  (47) 

is  of  great  interest.  It  is  a  nilpotent  2?-graded  Lie  algebra.  In  particular,  given  R  €  172^,  we 
can  use  equation  (44)  and  the  following  exact  sequence[31,  83] 

A  A^H.(R:Z!)  ^  ([iti.t.I  /  [».,[,r„x,]))  0  (48) 

where: 

xi  =  xi  (R) 

[  ,  ]  is  induced  by  the  commutator  and 
p.  is  dual  to  U:  A*  Hi  ►  H2  (R;^) 

to  obtain 


gr(xi)/xi*^  =  iri/[xi,xi]  =  Hi 


(49) 


«r(xi)  /  xp^  =  Hi  ©  Image([- ,  •])  (50) 

These  relations  have  already  been  considered  in  my  work  on  eddy  currents  in  multiply 
connected  regions.  In  order  to  consider  obstructions  to  Clebsch  potentials,  we  want  to 
consider 
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P(ir,(R))  /  (51) 

where  R  C  R  is  the  exterior  of  a  flux  tube  or  some  other  region  left  invariant  by  the  ‘‘flow’* 
of  curl  F.  The  literature  cited  in  the  bibliogTaphy[ll,  15,  16,  17,  18,  19,  20,  21,  22,  28,  29, 
33,  38,  39,  41,  42,  43,  58,  75,  77,  79,  80,  82,  83,  84]  makes  the  foUowing  connections: 


The  additional  algebraic  gadget  which  appears  in  three  dimensional  problems  is  Thurston’s 
norm  on  homology[87].  Its  relation  to  the  above  diagram  is  less  obvious  and  does  not  seem 
to  appear  in  the  literature. 


3.3  Clebsch  potentials  and  Hamiltonian  systems 

It  was  pointed  out  to  me  by  Professor  Marsden  [U.C.  Berkeley]  that  a  solenoidal  vector  field 
F  possessing  a  set  of  Clebsch  potentials  ^  and  17  implies  a  Hamiltonism  structure  for 


i  =  F(x) 

To  see  this,  consider  0(x)  and  compute 

^  ^  X  '^’i)  (“) 

or  ^  =  Vt;  •  (V^  x  V() 

Alternatively,  we  can  write 


(52) 


(53) 
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where  we  define 


(54) 


{a,b},  =  .  (Va  X  Vb)  (55) 

and  it  is  easily  verified  that  the  bracket  {  }  defines  a  Poisson  structure.  This  brii^  us 

back  to  the  theory  of  contact  structures.  That  is  the  global  theory  of  Hamiltonian  systems 
can  be  used  to  define  obstructions  to  finding  Clebsch  potentials. 

From  this  vantage,  it  is  clear  that  one  must  retreat  from  the  overly  ambitious  goal.  For 
the  antenna  synthesis  problem,  it  is  reasonable  to  look  for  additional  constraints  on  the 
current  distribution,  which  make  it  more  amenable  to  description  by  a  global  set  of  Clebsch 
potentials.  Unfortunately,  no  such  constraints  are  known.  Alternatively,  one  could  look  to 
other  models  of  physical  phenomena  involving  solenoidal  vector  fields  (e.g.  micromagnetics 
or  invisdd,  incompressible  fluid  flow)[l,  60]  to  find  simpler  Hamiltonian  systems  which  can 
guide  our  intuition. 
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4  Appendices 

4.1  Paper  published  in  which  APOSR  support  is  acknowledged 

“Metric  Dependent  Aspects  of  Inverse  Problems  and  Functionals  Based  on  Helidty” 

4.2  Previous  work  related  to  topological  accessibility 

“A  Topological  Invariant  for  the  Accessibility  Problem  of  Micromagnetics” 
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Metric  dependent  aspects  of  inverse  probiems  and  functionais 
based  on  heiicity 

P.  R.  Kotiuga 

Department  of  ECS  Engineering,  Boston  Uniuersity,  Boston,  Massachusetts  02215 

The  heiicity  of  a  vector  field  is  a  metric  independent  density.  Functionals  with  first  order  elliptic 
systems  for  Euler-Lagrange  equations  have  been  constructed  from  the  heiicity.  The  metric 
invariance  is  preserved  for  finite  element  discretizations  involving  “Whitney  elements.*’  This 
paper  relates  differential  geometric  aspects  of  inverse  problems  to  heiicity  based  functionals  in 
two  contexts.  First,  the  inverse  problem  of  electrical  impedance  tomography  in  isotropic  media 
is  known  to  be  equivalent  to  determining  a  metric  within  a  given  conformal  class  from  a  given 
“Dirichlet  to  Neumann”  map.  This  fact  is  related  to  the  heiicity  functional  and  Wexler’s 
algorithm  for  recovering  an  isotropic  conductivity.  Second,  Maxwell’s  equations  in  “spinor 
form”  are  shown  to  be  the  Euler-Lagrange  equations  of  some  complexified  time  dependent 
generalization  of  the  heiicity  functional.  In  this  case  metric  dependent  aspects  yield  insight  into 
the  “inverse  kinematic  problem  in  seismology.”  These  two  examples  illustrate  the  underlying 
geometric  structure  in  classes  of  inverse  problems  and  algorithms  for  their  solution. 

“Do  you  know  Grassmann’s  Ausdehnungslehre?  Spottiswood  spoke  of  it  in  Dublin  as 
something  above  and  beyond  quartemions.  I  have  not  seen  it,  but  Sir  William  Hamilton  of 
Edinburgh  used  to  say  that  the  greater  the  extension  the  smaller  the  intention.” 

“May  one  plough  with  an  ox  and  an  ass  together?  The  like  of  you  may  write  everything  and 
prove  everything  in  4nions,  but  in  the  transition  period  the  bilingual  method  may  help  to 
introduce  and  explain  the  more  perfect.” 

Excerpt  of  a  letter  from  J.  C.  Maxwell  to  P.  G.  Tait;  A.,  P.  Wills,  Vector  Analysis  with  an 
Introduction  to  Tensor  Analysis  (Prentice-Hall,  Englewood  Cliffs,  NJ,  1931),  pp.  XXV-XXVI. 


i.  INTRODUCTION 

The  issues  to  be  addressed  are  most  easily  introduced 
by  recalling  the  metric  and  coordinate  invariance  of  the 
heiicity  of  the  magnetic  field  intensity  vector  H.  Since  Am¬ 
pere’s  law  associates  the  magnetic  field  intensity  H  with 
integration  along  curves,  we  naturally  identify  it  with  a 
differential  1  form  as  follows: 

a=H-dt.  (1) 


Since 


■hji,  we  can  consider  the  A/y’s  to  be  the 


The  heiicity  is  then  defined  to  be 


ctiAdo) 


curl  H  dV, 


where  is  a  three-dimensional  region.  The  differential 
form  point  of  view'  is  inflicted  on  the  reader  at  the  onset 
smce  the  definition  (2)  is  clearly  metric  and  coordinate 
invariant.*  Furthermore,  these  properties  can  be  preserved 
in  the  context  of  the  finite  element  method  by  appealing  to 
Witney  forms.*  To  see  this,  consider  a  tetrahedron  A  with 
Vertices  denoted  by  pj,  l</<4,  and  define  edge  variables  h^j 
by 


jpi 


components  of  a  skew-symmetric  4x4  matrix  Given' 
these  edge  variables,  we  can  find  a  piecewise  linear  approx¬ 
imation  S  to  (U  by  appealing  to  Whitney  forms; 

4 

w=  Z  ffi^id^j.  (4) 

M=i 

Here  l</<4,  are  the  usual  barycentric  coordinates  on  A. 
When  discretized,  the  contribution  to  the  heiicity  from  one 
tertahedron  takes  on  the  amazingly  simple  form^: 

sj  aiAd«=Pfaff(^),  (5) 

where,  given  an  orientation  for  A,  Pfaff(<^)  is  a  square 
root  of 

Pfaff(<^)=Al2Aj4— A13/I24 4-^14^23  •  (fi) 

Equations  (5)  and  (6)  are  independent  of  interpolation 
node  coordinates  and  any  metric  dependent  bilinear  pair¬ 
ing  such  as  a  constitutive  law.  Furthermore  it  is  clear  that 
the  quadratic  form  associated  with  the  “element  stiffness 
matrix”  can  be  evaluated  with  only  three  multiplications! 
The  heiicity  functional* 

/(H,A)=  J^H(|curIH-l-/igradA-J)dF  (7) 
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,  ««plpits  the  metric  invariance  of  helicity.  Identifying  the 
current  density  J  with  a  differential  2  form,  the  helicity 
functional  is  seen  to  be  coordinate-independent  except  for 
the  term  involving  the  Lagrange  multiplier  A  which  en¬ 
sures  that  is  a  solenoidal  vector  field.  Again  this  prop¬ 
erty  is  preserved  when  the  functional  is  discretized  using 
Whitney  forms  and  the  discrete  Euler-Lagrange  equations 
have  a  remarkably  simple  interpretation  in  terms  of  Am¬ 
pere’s  Law.*  This  minimal  dependence  on  constitutive  laws 
and  coordinates  makes  the  Helicity  functional  ideal  for 
problems  requiring  iterative  solution  where  either  the  mesh 
formed  or  the  constitutive  law  is  updated  at  every 
ition.  This  aspect  will  be  developed  in  Sec.  II  when  we 
isider  Wexler’s  Algorithm  for  Impedence  Computed 
mography. 

To  introduce  a  complexified  time  dependent  helicity 
functional,  consider  first  the  Euler-Lagrange  equation  for 
the  functional  (7); 

/curl  H-fu  grad  A\  /J\ 

^(A)  =  (-div(^H/  )  =  (o)- 

If  we  choose  Cartesian  coordinates  and  set  /i=  1,  D  satis¬ 
fies 


(9) 

where  1  is  the  4x4  identity  matrix.^  Thus  D  is  an  elliptic 
partial  differential  operator  which,  in  some  concrete  sense, 
is  the  square  root  of  the  positive  definite  Laplace  operator. 
The  spectrum  of  D  is  the  square  root  of  that  of  the  Laplace 
operator  and  square  roots  are  distributed  in  a  manner 
which  is  nearly  symmetric  about  zero.  This  has  interesting 
consequences  for  the  numerical  solution  of  the  discretized 
equations.*  If  we  define  an  operator  D  by 

d 

D=D+iI-,  (10) 


then  using  Eq.  (9),  it  is  easy  to  verify  that 


DD*=I 


(11) 


Thus,  D  is  a  partial  differential  operator  which,  in  some 
concrete  sejse  is  the  square  root  of  the  wave  operator. 
Note  that  D  now  needs  to  operate  on  pairs  consisting  of  a 
complex  valued  vector  and  a  complex  valued  scalar.  Al¬ 
though  the  operator  D  is  closely  related  to  the  Dirac  op¬ 
erator  and  hence  the  quantum  mechanical  view  of  magne¬ 
tism,  our  purpose  is  to  extend  previous  work^  and  consider 
how  a  related  operator  reproduces  Maxwell’s  equations 
while  elucidating  the  coordinate  dependence  of  the  associ¬ 
ated  complexified  time  dependent  helicity  functional  whose 
Euler-Lagrange  operator  closely  resembles  D.  Thus,  Sec. 
Ill  of  this  paper  considers  a  complexified  time  dependent 
helicity  functional  whose  Euler-Lagrange  equations  are  re¬ 
lated  to  Maxwell’s  equations  just  as  Eqs.  (8)  are  related  to 
magnetostatics.  This  functional  is  then  considered  in  the 
context  of  inverse  problems. 


II.  WEXUER’S  ALGORITHM  REVISITED 

Consider  a  region  R  containing  a  given  inhomogeneous 
isotropic  conductivity  a.  Steady  current  flow  is  governed 
by 

divj=0,  J=aE,  curlE=0  in  /I.  (12) 

Letting  E  be  minus  the  gradient  of  the  scalar  potential  (^, 
we  can  prescribe  a  boundary  potential  on  dR  and  re¬ 
phrase  the  determination  of  J  in  terms  of  a  Dirichlet  prob¬ 
lem  for  (ft  : 

— div(agrad^)=0  in  R,  (13) 

(p=(^s  o"  (14) 

Given  a  solution  to  the  above  boundary  value  problem,  we 
can  calculate  the  Neumann  data  on  ■  grad  0  and  recover 
the  normal  component  of  J  on  dR.  Alternatively,  instead 
of  the  Dirichlet  data  (14),  J  can  be  uniquely  determined 
by  imposing  the  Neumann  data  on  dR.  From  a  variational 
point  of  view,  the  inhomogeneous  Neumann  data  is  diffi¬ 
cult  to  impose  and,  when  R  is  two-dimensional,  the  Neu¬ 
mann  problem  for  0  is  often  reformulated  as  a  Dirichlet 
problem  for  a  stream  function.  Similarly,  in  three  dimen¬ 
sions,  we  can  express  J  as  the  curl  of  a  vector  potential  T. 
Up  to  a  “gauge  transformation”  T  is  just  the  magnetic  field 
intensity  H  and  the  situation  is  analogous  to  the  use  of  a 
vector  potential  in  magnetostatics.  By  imposing  the  tan¬ 
gential  components  of  T  on  dR,  the  normal  component  of 
J  is  implicitly  prescribed.  When  J  and  E  are  eliminated  in 
Eq.  (12)  through  the  use  of  T,  current  boundary  condi¬ 
tions  are  imjx)sed,  and  the  “gauge”  of  T  is  fixed,  we  end  up 
with 

curl^-curlTj=0  in  R,  (15) 

div(aT)=0  in  R,  (16) 

nXT|=nXTa  on  dR.  (17) 

In  order  to  avoid  confusion,  we  call  Eqs.  (14)  and  (17), 
voltage  and  current  boundary  conditions,  respectively. 

For  a  given  inhomogeneous  conductivity  a,  there  is  a 
unique  correspondence: 

(ft  on  dR  mod  constants-*-*)?  ■  curl  T  on  dR  (18) 

which  we  call  the  voltage  to  current  map.  Impedence  to¬ 
mography  seeks  to  recover  a  given  this  voltage  to  current 
map.  In  practice,  given  a  collection  of  voltage  and  current 
measurements  {^gj,n  •  curlTg/}fLi,  we  would  like  to  re¬ 
cover  a  reasonable  estimate  for  <r  in  /?.  Reference  8  seems 
to  have  proposed  the  most  useful  algorithm  for  doing  this. 
See  Kohn  and  Vogelius’  for  a  survey  of  early  results,  and 
Sylvester  and  Uhlmann'°  for  a  geometric  understanding  of 
the  problem.  Our  concern  is  to  understand  how  the  helicity 
functional  points  to  a  useful  variant  of  Wexler’s  algorithm. 

Consider  the  helicity  functional; 

/(T,0)=  J^T  -  (j  curl  T-|-<r  grad  0)rfK  (19) 

The  Euler-Lagrange  equations  are 
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curl  T+a  grad  ^=0,  (20) 

_div(aT)=0.  (21) 

Taking  the  divergence  of  Eq.  (20)  yields  Eq.  (13)  while 
dividing  Eq.  (20)  by  a  and  taking  the  curl  yields  Eq.  ( 15). 
Equation  (21)  is  the  same  as  Eq.  (16).  Hence  any  pair 
(T,^)  satisfying  the  Euler-Lagrange  equations  of  the  he- 
licity  functional  (19)  automatically  satisfy  the  partial  dif¬ 
ferential  equations  (but  not  boundary  conditions!)  of 
steady  current  conduction  when  formulated  in  terms  of 
vector  or  scalar  potentials.  Furthermore,  recalling  Eq.  (9) 
we  see  in  what  sense  Eqs.  (20)  and  (21)  are  the  square 
root  of  the  combination  of  both  potential  formulations. 

We  now  demonstrate  why  the  metric  invariance  of  he- 
licity  is  so  useful  for  three-dimensional  impedance  tomog¬ 
raphy.  Suppose  we  have  the  data  •  curl  j  of  this 
inverse  problem  and  suppose  that  from  this  data  we  find  a 
consistent  set  of  pairs  of  boundary  conditions  for  the  he- 
licity  functional  {(^„,«  •  T„),(«XT„)}jIL,.  Furthermore, 
let  the  functional  (19)  with  principal  boundary  conditions 
(^„n  T,)  and  (nXT„)  be  denoted  by  and  I„i^,  re¬ 
spectively.  Given  an  initial  “guess”  Oq  for  the  conductivity, 
we  can  find  the  stationary  points  of  the  functionals 
and  where  l<«<Af.  In  general,  after  k  iterations  of 
Wexler’s  updating  scheme,  we  can  let  nnd 

b*  the  stationary  points  of  respectively. 

The  rule  for  updating  a  is  then  given  by 

XjL  1  curl T*,/  -grad 

grad  9k,y„  Vk,y^ 

Hence  every  iteration  of  the  proposed  algorithm  involves 
finding  the  stationary  values  of  IN  functionals 

l<n<.Ar  and  updating  a  by  the  rule  (22).  The  sur¬ 
prising  aspect  is  that  the  2N  functional  all  depends  on  an  in 
the  same  way!  Furthermore,  from  the  arguments  presented 
in  the  introduction  of  this  paper,  the  dependence  on  o  only 
occurs  in  the  coupling  of  ^  to  T  and  not  in  the  helicity  term 
of  the  functional.  For  finite  element  implementations,  this 
means  that  the  updated  a  given  by  Eq.  (22)  only  affects  a 
relatively  small  part  of  one  stiffness  matrix,  in  each  itera¬ 
tion! 


III.  HELICITY  AND  TIME  DEPENDENCE 

In  this  section  we  consider  a  time  dependent  version  of 
the  helicity  functional  and  the  relevance  of  metric  indepen¬ 
dence  in  the  solution  of  inverse  problems.  Let  y  be  an 
isotropic  constitutive  law,  (F,^)  a  pair  of  complex  func¬ 
tions,  one  vector  and  one  scalar,  and  consider  the  func¬ 
tional 


A:(F,^)=5(F,^)-f,!r*(F,^), 

where 


5(F,^) 


(23) 

dP) 


curlF  ydF  dP\ 

-+/2^+ygrad^-/^j 


d* 

+  — div  P  jy  d/.  (24) 


The  Euler-Lagrange  equation  for  A'(F,d)  is 


aF 

curl  F—iy— — \-y  grad  tf> 
at 

'  'a/  1 

—  div(yF) 

[-divPj 

(25) 


When  y=  1,  we  ^  that  the  operator  D  defined  by  Eq.  (25) 
is  the  operator  D  of  Eq.  ( 10)  if  it  were  not  for  a  missing 
time  derivative  of  (f).  This  is  an  intentional  omission  since, 
if  P  is  real,  Eq.  (25)  reduces  to  Maxwell’s  equations  when 


E+iH  =  F, 


(26) 


D-(-/B  =  yF, 


(27) 


(  ^P  \ 

(J,p)  =  ( -^,-divPj.  (28) 

With  these  identifications  the  conservation  of  charge  is 
satisfied  automatically,  the  variables  are  assumed  to  be  di¬ 
mensionless,  and  we  have  the  constitutive  laws  for  chiral 
media  if  y  is  complex.  When  ^  is  forced  to  be  zero,  Eqs. 
(25)  are  well  known.^’”  However,  for  nonzero  ^  there  are 
several  novel  aspects  of  the  present  formulation.  First,  ^  is 
a  Lagrange  multiplier  which  imposes  the  divergence  equa¬ 
tions^  That  is  Eq.  (25)  represents  eight  scalar  equations 
which  are  Maxwell’s  equations  in  component  form.  Sec¬ 
ond,  the  operator  D  represents  a  square  root  of  the  wave 
operator  applied  to  F  and  the  square  root  of  the  Laplace 
operator  applied  to  Third,  as  in  Sec.  II,  the  square  root 
operator  has  some  nice  metric  independent  aspects  which 
can  be  exploited  in  the  context  of  inverse  problems.  These 
follow  from  the  appearance  of  the  helicity  term.  The  in¬ 
verse  problem  connected  to  the  complexified  helicity  func¬ 
tional  is  the  inverse  kinematic  problem  of  seismology.  See 
Sec.  II  of  Ref.  10.  In  analogy  to  the  real  case,  the  metric 
independence  of  the  helicity  term  ensures  that  only  a  small 
part  of  the  functional  needs  to  be  updated  when  the  con¬ 
stitutive  law  is  updated. 

Support  from  AFOSR  is  gratefully  acknowledged. 
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A  topological  invariant  for  the  accessibility  problem  of  micromagnetics 

P.  R.  Kotiuga  and  R.  Giles 

Department  of  ECS  Engineering,  Boston  University,  Boston,  Massachusetts  022 15 

A  singularity-free  three-dimensional  micromagnetic  configuration  is  accessible  from  a  given 
initial  state  if  there  is  a  continuous  family  of  interpolating  states.  This  paper  proposes  an 
integral  formula  for  the  accessibility  invariant  which  classifies  (up  to  homotopy)  the  possible 
final  states  which  are  inaccessible  from  a  given  initial  state.  For  infinite  samples,  we  give 
explicit  formulas  describing  initial  states  from  which  the  final  state  of  uniform  saturation 
cannot  be  accessed. 


I.  INTRODUCTION 

Fundamental  to  explaining  the  operation  of  micromag¬ 
netic  memories  is  the  assumption  that  information  can  be 
stored  in  topological  configurations  of  magnetization  which 
can  be  made  energetically  stable.  The  notion  of  a  topological 
configuration  is  a  consequence  of  the  basic  tenet  of  micro¬ 
magnetics  which  asserts  that  the  magnetization  vector  field 
is  continuous  and  has  a  nonzero  length.  Unfortunately,  there 
are  circumstances  where  this  basic  tenet  must  be  relaxed  and 
one  considers  micromagnetic  singularities.  This  paper  intro¬ 
duces  the  notion  of  topological  accessibility  into  micromag¬ 
netics,  defines  an  accessibility  invariant  which  distinguishes 
mutually  inaccessible  states,  and  relates  the  accessibility 
problem  to  micromagnetic  singularities. 

The  previous  investigation  of  the  formal  properties  of 
Bloch  points'  plays  an  important  role  in  the  present  paper 
for  two  reasons.  First,  for  any  two  singularity-free  configura¬ 
tions  defined  at  times  f ,  and  ,  which  are  mutually  inacces¬ 
sible,  one  would  like  to  have  a  concrete  understanding  of 
how  and  why  Bloch  points  come  into  play.  In  the  case  where 
Bloch  points  cannot  enter  or  exit  through  the  boundary  of  a 
magnetic  sample  R,  Bloch  points  must  be  created  and  anni¬ 
hilated  in  pairs  of  opposite  “degree.”'  Hence,  it  is  possible  to 
regard  the  creation-evolution-annihilation  cycle  of  a  Bloch 
point  pair  as  a  closed  curve  in  i?  X  [  ,^2  ]  •  For  static  config¬ 
urations  this  helps  us  attribute  inaccessibility  to  reverse 
magnetized  tori  which  require  the  creation  and  annihilation 
of  a  Bloch  point  pair  in  order  to  collapse. 

The  second  reason  why  this  paper  builds  on  a  formal 
treatment  of  Bloch  points  is  that  the  “accessibility  invar¬ 
iant”  is  constructed  from  a  “vector  potential”  for  the  sole- 
noidal  gyrovector  field  and  the  motivation  for  this  construc¬ 
tion  is  most  transparent  through  the  use  of  differential 
forms.  Specifically,  Slonczewski’s  gyrovector  integral  for  the 
degree  of  a  mapping  into  the  Feldtkeller  sphere,^  when  re¬ 
phrased  in  terms  of  differential  forms,'  provides  the  lan¬ 
guage  which  relates  the  inaccessibility  invariant,  in  the  case 
of  infinite  samples,  to  the  Hopf  invariant.^  * 

Let  us  review  the  development  of  the  Hopf  invariant 
from  both  the  differential  form  and  gyrovector  points  of 
view.  Recall '  that  a  continuous,  nowhere  zero,  three-dimen¬ 
sional,  magnetization  vector  field  m  can  be  normalized  to 
have  unit  length  and  thought  of  as  a  mapping 

(1) 


where  R  is  the  magnetic  sample  and  5^  is  the  Feldtkeller 
sphere.  If  we  introduce  spherical  coordinates  (^,d)  on  5^ 
then  the  normalized  “volume  form”  on  (element  of  area 
divided  by  total  area)  is 

to=(l/47r)sinBtl0Ai/^.  (2) 

Definition  ( 1 )  makes  0  and  ^  functions  of  the  coordinates  on 
R.  We  can  “pull  back”  to  via  the  mapping/  to  obtain  a  differ¬ 
ential  2-form  on  R  in  terms  of  which  we  define  the  gyrovec¬ 
tor,  g: 

/*to  =  ( l/4tr)g-dS, 

g  =  grad(cos  0)  X grad  ( 3 ) 

Since  pull  backs  and  exterior  differentiation  commute 
(f*  ’d  sc  df* )  and  to  is  a.  closed  form  on  S  we  can  conclude 
that  the  gyrovector  is  solcnoidal.  That  is, 

dco  =  0=^0  =  d(f*(o)^divg  =  0.  (4) 

Assume,  for  the  time  being  that  the  magnetic  sample  R  occu¬ 
pies  all  of  R’  and  is  free  of  point  defects.  It  is  then  natural  to 
describe  the  gyrovector  by  a  vector  potential.  That  is,  we 
would  like  to  find  a  1-form  a  or  a  G  such  that 

da  —  f*(o,  i.e.,  curl  G  =  g.  ( 5 ) 

Assuming  that  m  is  sufficiently  close  to  the  constant  vector 
field  far  away  from  the  origin  in  R^,  the  integral 

/(/)=  f-  aAf*w  =  -~  \  GgdV  (6) 

Jr'  (4ir)^  Jr' 

is  finite  and  invariant  under  the  “gauge  transformation”: 
a-^a  +  dx,  i  e.,  G— G-l- grad;^',  (7) 

provided  that  x  is  close  to  constant  far  away  from  the  origin. 

The  integral  defined  by  Eq.  (6)  always  takes  on  integer 
values.  It  is  the  purpose  of  this  paper  to  show  how  the  var¬ 
ious  nonzero  integral  values  of  (6)  classify  the  micromag¬ 
netic  configurations  which  are  inaccessible  from  the  field 
consisting  of  a  constant  m  throughout  R^  and  to  generalize 
the  construction  to  finite  samples  R  and  nonconstant  initial 
or  final  states.  Traditionally,  the  Hopf  invariant  is  defined 
for  mappings  of  S '  ( the  unit  sphere  in  R^ )  onto  S  but,  as 

we  shall  see,  R^  is  easily  mapped  onto  by  stereographic 
projection  and  the  asymptotic  conditions  on  m  (and  hence, 
g)  ensure  that / can  be  regarded  as  a  map  from  5 '  to  S 
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<3  be 

FIG.  1.  A  reverse  magnetized  torus. 


II.  THE  HOPE  INVARIANT  AND  MICROMAGNETICS 


Given  a  magnetization  configuration  m(r,r)  on  all  of 
at  some  initial  ,  we  can  ask  whether  there  is  a  singularity- 
free  time  evolution  such  that  at  time  ,r)  is  some  state 

of  uniform  saturation  given  by  a  constant  vector  m^.  For 
example,  consider  the  reverse  magnetized  torus  of  Fig.  1(a) 
where,  in  Feldtkeller’s  notation,’  Fig.  1(c)  represents  the 
cross  section  depicted  in  Fig.  1  (b).  A  simple  rescaling  argu¬ 
ment  shows  that  shrinking  the  torus  makes  m  (/,r )  discontin¬ 
uous  in  time  at  (2  ■  Assuming  that  the  torus  would  like  to 
shorten  its  length  without  changing  its  overall  diameter,  we 
expect  it  to  “break  symmetry”  and  collapse  by  creating  a 
Bloch  point  pair  which  is  annihilated  when  the  final  state  is 
achieved.  In  general,  the  integral  of  Eq.  (6)  gives  an  integer 
whose  absolute  value  is,  in  some  vague  sense,  the  minimum 
number  of  Bloch  point  pairs  required  to  be  created  and  anni¬ 
hilated  in  order  to  obtain  the  state  of  saturation.  In  the  case 
of  Fig.  1  we  shall  see  that 

16ir*  Jr' 


In  order  to  make  this  plausible,  we  need  some  formal  results. 
Consider  the  unit  sphere  5’  in  R'*  to  be  the  unit  sphere  in 
defined  by  2,z,  =  1.  Stereographic  projection,  de¬ 

noted  by  cTj ,  then  maps  R’  onto  S  ’  where 


<7,:  (X|,JCj,a:3)-(7,,Zj) 

_  [2jC|  -I-  21x2,2x2  -f  /( 1  -  7:7  -  Xj  -  Xj )  ] 
1  -I-  xf  4-  X j  -t-  xj 

(7,“':  (Z,,7j)^(X,,Xj,X3) 


(9) 


(Rez,  ,Imz,  ,ReZ2) 
1  —  Im  Z2 


when 


Im  Z2 


(10) 


Hence,  when  m  approaches  the  constant  field  as 
x^t  +  xl  +  xl 00  there  is  no  loss  of  generality  in  consider¬ 
ing  the  magnetization  vector  as  a  mapy£73"  '  from  S’ to  5^.  In 
this  case,  Eq.  (6)  is  the  traditional  Hopf  invariant.  The  main 
results  concerning  the  Hopf  invariant  are  summarized  by  the 
following.’'^ 

Theorem  1:  (a)  The  Hopf  invariant  is  independent  of 
the  choice  of  a  (i.e.,  G)  in  Eq.  (5).  (b)  Homotopic  maps  / 
have  the  same  Hopf  invariant,  (c)  If  two  maps,y573~  ‘  and 
faf  '  have  the  same  Hopf  invariant  then  they  are  homotopic 
in  S’. 

Here  homotopic  means  that  there  is  a  continuous  family 
of  interpolating  states.  For  the  above  example  of  a  reverse 
magnetized  torus,  the  configuration  is  homotopic  to  the 
Hopf  fibration  for  which  a  straightforward  calculation’ 
shows  that  I{J^ )  =  1.  That  is,  Eq.  (8)  is  indeed  correct. 

Frequently,  the  integral  expression  given  by  Eq.  (6)  is 
too  cumtersome  to  be  an  effective  computational  tool.  One 


then  falls  back  on  the  linking  number  interpretation  of  the 
Hopf  invariant.’  Consider  any  two  regular  values p,  and  p2 
for  the  map /.  That  is,  the  Jacobian  of /  has  rank  2  at  p^  and 
P2  ■  The  preiraages  of faf  ‘  of  these  two  points  are  “compact 
manifolds”  and,  hence,  disjoint  circles  in  5  ’.  The  Hopf  in¬ 
variant  is  the  linking  number  of  these  two  circles: 

/(/)=;  Link[a/-'(p,), <7/- '(Pj)].  (H) 

The  proof  of  the  equivalence  of  Eqs.  (6)  and  (11)  makes 
clear  why  the  Hopf  invariant  is  always  an  integer.  Returning 
to  Fig.  1  we  see  that  Link[/r  '(  —  )/,”  '(  4-  )  ]  is  equal 

to  one.  Hence,  the  reverse  magnetized  torus  cannot  evolve 
without  singularities  into  the  state  of  uniform  saturation. 

We  now  consider  a  more  explicit  treatment  displaying 
states  with  higher  Hopf  invariant.  A  family  of  maps  ,  from 
R’  to  S  ’  describing  a  continuous  nowhere  vanishing  magne¬ 
tization  field  m  is  given  as  a  composition:^  ‘ 

4,  =  r^2Ag,,<73,  (12) 

where,  introducing  coordinates,  we  write 
R’--5’-5’-CP‘-.S^ 

*»,  h  ”1 

(x,,Xj,Xj)-(z,,Z2)-(u;,,u)2)-u)  -(m,  4-/m^,m,). 

Here  S’  is  the  unit  sphere  in  C’{fP*  is  the  complex  projective 
space,  CTj  and  czj  are  stereographic  projection  into  three-  and 
two-dimensional  spheres,  respectively,  and  aj  h  is  the  Hopf 
fibration.  Explicitly,  is  given  by  Eq.  (9),  while  the  re¬ 
maining  functions  are  defined  as  follows: 


Sp., :  (Zi  ^2 )  -'  (z?  ^ )  =  ( u'l  .U'2 )  with  p,9eZ; 
A:  (lu,  ,u>2 )  -  /Wi  —  w; 


w- 


{2w,l  -  |u;|^) 


=  (m„  4- 


(13) 


1  4-  1^1^ 

When  p  =  1  =  9  the  composite  map  is  homotopic  to  the  con¬ 
figuration  described  in  Fig.  1  and  from  the  discussion  of  Sec. 
II,  we  conclude  that  /((72  A)  =  1.  Since  the  Hopf  invariant 
can  be  described  in  terms  of  the  degree  of  a  map,’  it  follows 
that 


/(/^CTj  ' )  =  )  =  deg(gp., ).  ( 14) 


Also,  since  the  function  zT  has  degree  m  when  regarded  as  a 
map  from  the  unit  circle  in  the  complex  plane  to  itself,  it 
follows  from  Eq.  ( 14)  that 


')=^^g(gp.p)=P(l-  (15) 

Thus  from  theorem  1  we  conclude  the  following: 

Theorem  2:  (a)  Two  configurations^.,  defined 

by  Eq.  (12)  are  mutually  inaccessible  if  pq^p'q'-  (b)  If 
pq  =  p'q'  then  the  maps  and  are  mutually  accessible 
in  S’  but  the  map^.,  may  not  be  accessible  from 4,,.  in  R’ 
using  a  finite  amount  of  energy. 

To  obtain  a  concrete  and  intuitive  appreciation  of  these 
results,  it  is  useful  to  obtain  explicit  formulas  for  the  maps 
4.,  in  terms  of  torroidal  coordinates  on  R’  (Ref.  7): 


(X,  +  1X2  ,Xj) 


_ 1 _ 

coshfu)  —  cos(t/) 


X  [sinh(i/)e'*,sin(M)l,  (16) 
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where  0<^<2jr,  0<,u<.2ir.  Substituting  Eq.  (16) 

into  (9)  and  simplifying  we  obtain 

(^1^2)= — -  (sinh(i')e'V“]-  (17) 

cosh(i/) 

Equations  (13)  then  give 


(m,  + 


(Zae-*/-* -■>“>,  1  -g^) 

1+a^  ’ 


a  =  sinh^(i/)cosh‘'  '’(f). 


(18) 


Equation  (18)  describes,  in  terms  of  torroidal  coordinates, 
continuous,  nowhere  zero  tr"gnetization  distributions  with 
Hopf  invariant  pq.  The  special  case  ofp  =  q  =  1  is  described 
in  Fig.  1.  In  general  there  is  a  lirp  twist  in  the  direction  of 
increasing  ^  and  a  lirq  twist  in  the  direction  of  increasing  u. 
Returning  to  the  linking  number  interpretation  of  the  Hopf 
invariant,  the  generic  curves  of  constant  m,  with  linking 
number  pa,  are  found  by  fixing  p(f>  —  qu  and  v. 


III.  THE  ACCESSIBILITY  INVARIANT 


We  now  generalize  the  Hopf  invariant  to  obtain  a  topo¬ 
logical  accessibility  invariant  which  is  useful  in  the  context 
of  finite  samples  for  which  the  magnetization  is  given 
throughout  the  sample  at  initial  and  final  times  and  on  the 
boundary  of  the  sample  as  a  function  of  time.  Consider  a 
finite  region  R  and  a  time  interval  [  r,  ,t2  ]  denoted  by  I.  The 
region  R  is  assumed  to  have  no  defects  or  “holes’*  so  that  the 
( CO )  homology  groups  are  trivial.  The  magnetization  on  the 
boundary  d(R  X/)  of  the  region  of  space-time  is  fixed  by 
specifying  m(f,  ,r)  and  m(t2,r)  for  ail  reR  and  m(/,r)  for 
t^  </<r2  and  re  d  R.  Here  it  is  natural  to  regard  time  as  the 
parameter  describing  the  homotopy  between  mutually 
accessible  initial  and  final  states.  The  situation  is  as  before 
except  that  a  continuous,  nowhere  zero  magnetization  field, 
if  it  exists,  would  be  described  by  a  map  f  which  extends /: 

f-.R  Xl-~S\  where  /:d(R  X/)  ( 19) 

If  R  has  trivial  (co)homology  groups  then  it  is  possible  to 
find  a  potential  1 -forms  a  and  5  to  describe  the  pull  back  of 
the  normalized  volume  form  on  5’: 


da=f*aj  on  R  X/,  da=f*co  on  d(R  X/),  (20) 
Consider  the  integral  defining  the  accessibility  invariant: 

J{J)  =  j  a/\da.  (21) 

Jam  'It 

The  value  of  this  integral  is  zero  if  the  map / can  be  extended 
to  a  singularity-free  map  /  because  if  such  an  /  exists,  by 
Stokes’  theorem  and  Eq.  (20)  we  have 

I  df\da={  d(aAda) 

Ji>(  R  y  l)  Jr  </ 


=  r  da  Add  =  f  (7*(«u)  A/*(<u)  ] 
Jr  y  I  Jr  y  I 


(a) Aw)  =0 since <y Aw  =  0 on 


Hence  a  nonzero  value  forJif)  expresses  the  impossibility  of 
extending  the  data  given  on  d(.R  X I)  to  a  continuous  non¬ 
zero  distribution  of  magnetization  on  R  X/.  As  in  the  case  of 
the  Hopf  invariant,  Stokes’  theorem  is  used  to  verify  that  the 


integral  J(f)  is  invariant  under  the  gauge  transformation 
a — a  -H  dx-  The  linking  number  interpretation  is  also  valid 
so  that  the  integral  /(/)  takes  on  integral  values  which  repre¬ 
sent  the  linking  number  of  two  preimages  of /in  the  closed 
manifold  d(R  X/).  Indeed  we  expect  such  a  result  since  Eq. 
(21 )  reduces  to  Eq.  (6),  when  m  is  a  constant  vector  on  dR 
for  all  time,  over  R  at  .  and  on  R’  —  R  at  r, . 

A  concrete  understanding  of  Eq.  (2 1 )  in  terms  of  vector 
and  scalar  functions  on  X  /  is  obtained  by  appealing  to  an 
analogy  with  electrodynamics.  The  electromagnetic  field  is 
described  by  a  differential  2-fonn  P  on  space-time  which 
being  closed  can  be  expressed  as  the  exterior  deri  .  ative  of  a 
potential  1-forma: 

da=P  =  BdS  +  Edldt, 
where 

a  =  Adi-<j>dt.  (22) 

In  micromagnetics,  P  plays  the  role  of f*w  so  that  the 
vector  fields  A  and  B  are  formally  analogous  tc  G  and  g. 
Thus,  to  complete  the  prescription  given  by  Eq.  (20)  we 
need  to  introduce  an  analog  of  ^  and  E.  Introducing  a  scalar 
potential  y,  the  equation  da  —  P  shows  that  Eq.  (20)  can  be 
rewritten  as 

da  =  curl  G  dS  +  ^grad  y  —  ^j^  dlAdt  =  f*w.  (23) 


Once  a  potential  (G,y)  is  found  Eq.  (21 )  becomes 

^  L  H  r  -  f  I  -  l-  curl  -i' 

-fj  G  curlCItj  dK- J  G  curlGl/, 


Invariance  of  Eqs.  (23)  and  (24)  under  gauge  transforma¬ 
tions  of  the  form 


G-C-f  gradjt',  y--y  +  ~  (25) 

dt 

are  still  best  understood  in  the  context  of  equation 
“a  — a  -f  dx"-  The  results  of  theorem  1  and  the  calculation 
below  Eq.  (21)  are  stated  as  follows: 

Theorem  3:  The  integral  given  in  (21 )  or  (24)  takes  on 
integer  values  which  are  nonzero  only  if  the  initial  and  final 
states  given  in  ( 19)  are  mutually  inaccessible.  If  R  is  a  con- 
tractiole  region  of  finite  extent,  the  states  are  topologically 
accessible  whenever  the  integral  has  a  zero  value. 

The  relevance  of  Eq.  (6)  to  micromagnetics  has  also 
been  noted  by  P.  Asselin,  who  provided  the  example  illus¬ 
trated  in  Fig.  1.  The  authors  would  like  to  thank  Honeywell 
Inc.  for  financial  support  of  this  work  and  Floyd  Humphrey 
for  his  encouragement. 
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